Density functional theory study of (OCS)2 



G. Bilalbegovic 

Department of Physics, Faculty of Science, University of Zagreb, Bijenicka 32, 

10000 Zagreb, Croatia 



Abstract 

The structural and electronic properties of the carbonyl sulfide dimer anion are 
calculated using density functional theory within a pseudopotential method. Three 
geometries are optimized and investigated: C2v and C2 symmetric, as well as one 
asymmetric structure. A distribution of an excess charge in three isomers are stud- 
ied by the Hirshfeld method. In an asymmetric (OCS)^ isomer the charge is not 
equally divided between the two moieties, but it is distributed as OCS"'^'^- OCS"'^'^. 
Low-lying excitation levels of three isomers are compared using the time-dependent 
density functional theory in the Casida approach. 



1 Introduction 

Ionic nanoparticles of various sizes offer possibilities to investigate the change 
of a chemical bonding from molecules and clusters to the bulk. Studies of ionic 
clusters are important for understanding of microscopic aspects of chemical 
reactions in controlled environments [1]. Carbonyl sulfide may play an impor- 
tant role in an explanation of the origin of life on Earth. It acts as a condensing 
agent in the formation of peptides from amino acids [2]. This material is known 
as a component of volcanic gases and it is present in the atmosphere of Earth 
where it is a part of the global sulfur cycle. Carbonyl sulfide is also present in 
the atmosphere of Venus and in many other astrophysical objects. Therefore, 
OCS and related ions and clusters are important for various chemical and 
physical processes. It has been found that OCS~ anion is metastable, but that 
an addition of OCS or H2O molecules produces a stable structure [SP] . 

For homogeneous anionic clusters, such as (OCS)~, it is important to find out 
whether an excess electron is located in a vicinity of a particular monomer, or 
it is shared between two or more fragments. The photodissociation dynamics 
of (0CS)2 has been studied in the ion beam apparatus using a tandem time-of- 
fiight mass spectrometry [5] . The photoelectron imaging method has been also 
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applied in the studies of (0CS)2 Pl6] . The coexistence of two covalent OCS2 
isomers (having C2,; and C2 symmetry), as well as electrostatically bound 
OCS~- OCS have been proposed [Tlf5] . The formation of (OCS)~ cluster ions 
have been studied using low-energy electron attachment to molecular clusters 
(OCS)„ seeded in helium carrier gas [7j. This experiment have been supported 
by the coupled cluster theory calculations of electron affinity and vertical 
detachment energy of the carbonyl sulfide [7]. 

Advanced experimental methods, such as photoelectron spectroscopy, are able 
to explore various nanostructures. It is important to compare these results 
with different theoretical simulations of the same systems. The properties of 
(0CS)2 cluster have been studied using the Gaussian quantum chemistry code 
in connection with the photodissociation experiments A related (OCS)^ 
complex has also been investigated using the Gaussian program [8] . The pseu- 
dopotential based density functional theory is characterized by the computa- 
tional efficiency and recent intensive developments. It is nowadays possible to 
use this method for big systems which exhibit complex chemical and physical 
processes. It allows computations of larger systems than traditional quantum 
chemistry methods. The pseudopotential time-dependent density functional 
theory is promising and still in development for many applications, for ex- 
ample in studies of charge transfer processes. It is important to use density 
functional theory and pseudopotentials to simulate some of the structures 
and processes studied in experiments on time-resolved dynamics and solvation 
in big anion clusters [1]. However, methods initially developed in theoretical 
chemistry sometimes better describe dynamics of electrons involved in chem- 
ical processes. Therefore, it is important to test the pseudopotential density 
functional methods in studies on small charged clusters. In the previous com- 
putational study of (OC5')2 , which has been done using the Gaussian package 
[T][5] . the charge population has not been analyzed, and a direct calculation of 
low-lying excited energy levels of different isomers has not been presented. In 
this work (OCS)^ is studied using density functional theory (DFT) electronic 
structure calculations. Three geometries of the carbonyl sulfide dimer anion 
are minimized and their structures, energies, the ground state electronic prop- 
erties, as well as distributions of the lowest excitation levels are investigated. 



2 Computational Method 

DFT based calculations have been recently used as a powerful method for stud- 
ies of bulk, surfaces, and nanostructures of various materials [H]. The structural 
and electronic properties of (0CS)2 are studied in this work using the plane 
wave density functional Abinit code [10]. This code is robust, accurate and 
widely used. It has been applied by many users on various geometries and 
processes, and it is therefore well tested. The local density approximation of 
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the exchange functional in the parametrization of Perdew and Wang is apphed 
[TT] . The (0CS)2 ions (and related structures used for testing and compar- 
ison) are positioned in the middle of the periodically replicated simulation 
box. The side of the box for most simulations was 30 a.u. The Brillouin zone 
is sampled using the F point. Relaxation of ions is carried out by performing a 
series of self-consistent calculations. The Broyden-Fletcher-Goldfarb-Shanno 
minimization method is used for a structural optimization and all atoms are 
allowed to move. 

In the Abinit package charged systems are immersed in a neutralizing jellium 
background. For such charged particles artificial electrostatic interactions be- 
tween the initial system and its periodic images in the supercells are possible 
p!2|13fl4j . However, the supercells used here are sufficiently large to avoid 
these effects. To check this point several calculations for larger cells (35 and 
40 a.u.) are done and the differences in the total energy, bond lengths and 
angles are found to be 0.1% at the most. In addition, the charge density plots 
do not change when the side of the cell increases. 

Because of a fragility of some possible (0CS)2 structures, two kinds of pseu- 
dopotentials are used and the results are compared. First, the calculations are 
carried out using the TrouUier-Martins (TM) pseudopotentials [15], prepared 
by the Fritz Haber Institute code [16] and taken from the Abinit web page 
[To] . The cutoff of 35 Ha is used for the TM pseudopotentials. Calculations 
are also done using the relativistic separable pseudopotentials in the form of 
Hartwigsen, Goedecker and Hutter (HGH) [T7]. The cutoff of 50 Ha is applied 
for the HGH pseudopotentials. It has been found that the properties of CS, 
CS2 and several oxide molecules calculated using these pseudopotentials agree 
with experiments pT]. It is calculated in this work that both kinds of pseu- 
dopotentials (HGH and TM) produce a good agreement with experiments for 
the properties of OCS. This molecule (in experiments, as well as in calculations 
described here) stabilizes into a linear structure. The S-C distance is 1.548 A 
(the experimental value is 1.561 A P]), and the C-0 distance is 1.160 A (the 
experimental result is 1.156 A) in the calculations using the HGH pseudopo- 
tentials. The calculated distances are 1.548 A and 1.148 A, respectively for 
the S-C and C-0 bonds, when the TM pseudopotentials are applied. For the 
(0CS)2 isomers, only small numerical differences are found between calcula- 
tions carried out using the TM and HGH pseudopotentials. It is important 
to report the result for an electron affinity of the OCS molecule. However, it 
is not possible to stabilize (OCS)~ structure within the method used in this 
work and this is a necessary step for the calculation of electron affinity. The 
absence of a stable (OCS)~ anion in this calculation is in agreement with the 
experimental results. This anion is not formed in a standard ion source and 
almost no (OCS)~ is detected in the mass-spectrum [T|3llH[5|6] . 

The excitation energies are calculated within the time-dependent local density 
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Fig. 1. The density functional optimized structures and distances of symmetric 
(0CS)2 : (a) C2d, (b) 02- Bond lengths (in Angstroms) and angles (in degrees) are 
shown. 

functional theory using the Casida approach p^p9] implemented in the Abinit 
code. In this method the time- dependent density functional theory (TDDFT) 
equations are studied in the frequency domain. The eigenvalue problem is 
solved for the matrix whose main part is formed of the squares of differences 
between occupied and unoccupied Kohn-Sham electronic energies. The cou- 
pling matrix which includes the contribution of the Coulomb and exchange- 
correlation interactions is also added to the Casida matrix. It is known that 
low-lying excited states of small, chemically simple clusters and molecules are 
well represented in the Casida approach. New algorithms in TDDFT are neces- 
sary for periodic systems and high-lying excited levels. In this work the Casida 
approach to TDDFT is used to compare low-lying excitation energy levels of 
three (0CS)2" isomers. 



3 Results and discussion 

Three equilibrium structures of {0CS)2 are stabilized within approximations 
of density functional theory used in this work and for both kinds of pseudopo- 
tentials. These calculations started from initial geometries close to those pro- 
posed by Sanov and coworkers pl5| . Several other initial structures of (0CS)2 
were also considered, but none of them is found to be stable. The equilib- 
rium structures of (OCS)^, as well as corresponding bond lengths and angles 
calculated using the TM pseudopotentials are shown in Figs. 1 and 2. In the 
asymmetric (OCS)^ structure (shown in Fig. 2) both OCS fragments are non- 
linear. This also applies for the same isomer optimized using the HGH pseu- 
dopotentials. The electron charge isosurface plot for the asymmetric (0CS)2 
structure is also shown in Fig. 2. Two moieties are bonded between the carbon 
atom of one fragment and the sulphur atom of the other. The distance be- 
tween these two atoms is 2.49 A. Figures 2 and 3(c) show that a week covalent 
interaction between the two moieties exists. In the structures found by Sanov 
and coworkers one isomer is asymmetric. This is an electrostatically bound 
structure where one of the moieties is linear, and the charges are distributed 
as OCS~- OCS [TlfS] . Both asymmetric isomers, weakly covalent and electro- 
statically bonded one, may form under experimental conditions. The binding 
energies are presented in Table 1. The C2v structure is the most stable. The 
differences in energy of about AE'i = 0.25 eV and Ai?2 = 0.61 eV exist be- 
tween the C2,; structure and two other isomers. These results also show that 
the binding energy of the most stable structure of the C2^ (OCS)^ isomer is 
~ 3 eV lower than one of the neutral G2V (0CS)2 cluster optimized by the 
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Fig. 2. The density functional optimized asymmetric (0CS)2 structure. The electron 
charge density, as well as bond lengths (in Angstroms) and angles (in degrees) are 
shown. The charge contour corresponds to a constant density equal to 0.02 eV/ag. 

Table 1 

Binding energies (in eV) of three optimized anions of (0CS)2^ and a neutral (0CS)2 
cluster (having C2v symmetry) calculated using the TM and HGH pseudopotentials. 



Structure 


C2. (0CS)2 


C2 (0CS)2 


asymmetrical (0CS)2 


(0CS)2 


TM 


-46.794 


-46.544 


-46.188 


-43.750 


HGH 


-46.602 


-46.357 


-46.002 


-43.293 



Fig. 3. The HOMO states of (0CS)2 isomers: (a) symmetric, (b) C2 symmetric, 
(c) asymmetric structure. The light (dark) isosurface corresponds to the negative 
(positive) part of the wave functions. 

same method. 

The Hirshfeld method is a suitable technique to analyze the charge density 
distribution in molecules decomposed into atomic fragments [201121] . In this 
method a molecule is divided into atomic contributions proportional to the 
weight of the free atom charge density. The free atoms situated in their cor- 
responding positions in the molecule define the promolecule. This reference 
state of the promolecule is compared to the actual charge on the real molecule 
induced by the formation of chemical bonds. The Hirshfeld charges of (0CS)2 
are shown in Tables 2 and 3. The third structure (the results shown in Table 
3) is not symmetric and charges are different for two atoms in three pairs 
of the same kind (S, C, and O). A small positive charge is located on the 
carbon atoms in C2V and C2 symmetric structures, whereas a negative charge 
is distributed on all atoms in asymmetric (0CS)2 • The results for the neu- 
tral (0CS)2 dimer (also presented in Table 2) show that charges on all atoms 
become more negative in (008)2- The oxygen Hirshfeld charges in (00S)2 
isomers are substantially lower than in the water molecule and in the carbonyl 
oxygen of the acetamide, but it is also higher than in the carbonyl oxygen of 
the formaldehyde [22]. These Hirshfeld charges for various molecules are cal- 
culated within different ab initio computational programs. However, previous 
calculations (for a different material) have been shown that the values of Hir- 
shfeld charges are not sensitive to the calculation technique |[23j . Table 3 shows 
that the excess charge in an asymmetric (OOS)^ structure is distributed as 
—0.6 : —0.4. More bended fragment of this isomer takes more, but not all 
of the negative charge. The part of the excess charge is located on the other 
moiety which is also bended. 

Figure 3 presents the HOMO wave functions, whereas Fig. 4 shows the energy 
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Table 2 

Hirshfeld charges of (0CS)2 for two symmetric structures calculated using the TM 
pseudopotentials. The charges on atoms in the C2v symmetric neutral (0CS)2 clus- 
ter are also shown. The values for the HGH pseudopotentials are in the parenthesis. 



Structure 


C2. (ocs)2- 


C2 (0CS)2 


C2. (0CS)2 


6Q{S) 


-0.251 (-0.250) 


-0.296 (-0.294) 


+0.049 (+0.049) 


6Q{C) 


+0.026 (+0.018) 


+0.038 (+0.030) 


+0.100 (+0.092) 


SQiO) 


-0.274 (-0.268) 


-0.241 (-0.235) 


-0.149 (-0.141) 



Table 3 



Hirshfeld charges for the asymmetric (OCS)^ structure calculated using the TM 
pseudopotentials (the values for the HGH pseudopotentials are in the parenthesis). 
Atoms labeled with "1" belong to the moiety on the right in Fig. 2. 



Structure 


Atom 1 


Atom 2 


6Q{S) 


-0.156 (-0.154) 


-0.344 (-0.343) 


6Q{C) 


-0.006 (-0.016) 


-0.010 (-0.016) 


SQiO) 


-0.232 (-0.226) 


-0.245 (-0.238) 



levels of optimized (0CS)2" isomers. Delocalized negative orbitals exist only in 
the C2 symmetric (Fig. 3 (b)) and asymmetric structure (Fig. 3(c)). The cr-like 
bonding orbital exists between the carbon atom on the left fragment and the 
sulphur atom on the right one in the asymmetric structure (Fig. 3(c)). The 
details of electronic structure, such as a distribution of levels and shapes of 
HOMO states, are different for these three isomers. As in many other systems, 
these results show that the chemical and physical properties of nanostructures 
are determined not only by their size and the atoms they consist of. The details 
of geometrical and electronic structure are also important. 

In experimental studies of (008)2" [3ll5|6] it has been found that covalent iso- 
mers and asymmetric electrostatically bound OCS~- OCS structure behave 
differently in photoexcitation. The availability of low-lying excited states has 
been suggested to explain the autodetachment electron emission from the co- 
valently bound isomers. The results in Fig. 5 show that the most pronounced 
low excited levels exist for the C2 symmetric isomers. The low-lying excited 
states do not appear in the asymmetric (0CS)2 structure. Therefore, this cal- 
culations shows that the C2 covalently bound isomers are the most suitable 
for the autodetachment mechanism of electron emission. 

In summary, three electronic isomers of the carbonyl sulfide dimer anion are 
stabilized and investigated using the pseudopotential based DFT computa- 
tional methods. These three isomers of (OCS)^ have different geometrical 
and electronic properties, and therefore their participation in various chem- 
ical and physical processes is different. In particular, the Hirshfeld method 
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Fig. 4. Electron eigenvalues for (0CS)2 isomers: (a) C2v symmetric, (b) C2 symmet- 
ric, (c) asymmetric structure. Several levels are degenerate and almost degenerate, 
and therefore they are not distinguishable in the figure. 

calculations show that a distribution of a negative charge strongly depends on 
the structure and that in an asymmetric isomer the charge is distributed as 
—0.6 : —0.4, where the more bended moiety takes more charge. The low-lying 
excited levels of (0CS)2 isomers are calculated within the time-dependent 
density functional method and the lowest excitations are found for the C2 
symmetric structure. This work describes three possible structures of the car- 
bonyl sulfide dimer anion. These results are in rather good agreement with the 
previous computational results obtained by different method p!|5] . Two similar 
covalent isomers are obtained, but electrostatically bound OCS~- OCS struc- 
ture is not found in the present study. Instead weak covalent OCS~°'^- OCS~°'^ 
asymmetric structure is minimized. Many isomers, stable and metastable ones, 
may form under experimental conditions. Therefore, the existence of all four 
isomers are possible in experiments, and all of them may contribute to the 
time-resolved dynamics and solvation processes in {OCS)~ clusters. It is pos- 
sible to assess the reliability of these calculations by comparison with the 
results of the Gaussian code [T|5] . The two lowest energy covalent isomers of 
the same symmetry are found using both methods and their computationally 
optimized structural parameters are close; the largest differences are below 
2%. This agreement shows the precision and the reliability of both techniques, 
as well as the strengths of all modern computational methods in the studies 
of materials. The details of the structure and charge distribution in isomers 
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Fig. 5. The lowest excitation energies from the ground-state of (0CS)2 isomers: (a) 
C2v symmetric, (b) C2 symmetric, (c) asymmetric structure. The first ten mixed 
singlet and triplet energy levels are shown (several levels are not distinguishable on 
these figures). 

of (OCS)^ are also important on a more general level because they play a 
role as a prototype for the (OCS)~ cluster ions and corresponding solvation 
effects. Because of the computational efficiency of the pseudopotential den- 
sity functional theory method, it is possible to extend these calculations to 
larger anion clusters, such as {OCS)^ already studied in experiments described 
in Refs. p1l3]|5] . The pseudopotential time- dependent density functional tech- 
niques, used in this work to compare distributions od low-lying excited energy 
levels, are now under intensive development. These methods should be useful 
in an interpretation of experimental results on the time-resolved dynamics of 
cluster anions and other nanostructures. 
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